Heterogeneity in Classical and Non-Classical Nucleation 
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We investigate heterogeneous and homogeneous nucleation in nearest-neighbor and long-range 
Ising models for various quench depths. We find that the system has a true crossover from hetero- 
geneous to homogeneous nucleation for increasing quench depth only if the interaction is sufficiently 
long-range. The survival curves, defined as the fraction of systems that remain in the metastable 
state after a given time, have qualitatively different shapes as a function of quench depth for hetero- 
geneous and homogeneous nucleation when the interaction is short-range, but have identical shapes 
within the accuracy of our data for long-range interactions. 



Nucleation, the process by which a metastable state 
decays, plays an essential role in a wide variety of sys- 
tems. The nucleation mechanism involves the appearance 
of a critical droplet which overcomes a free energy barrier 
and initiates a decay into the stable state. In homoge- 
neous nucleation the droplet forms due to spontaneous 
fluctuations. Heterogeneous nucleation occurs when the 
droplet forms with the help of a wall, defect, or an im- 
purity such as an aerosol. Both homogeneous and het- 
erogeneous nucleation are technologically important, but 
much greater progress has been made in understanding 
homogeneous nucleation Q, H, IE Q than heterogeneous 
nucleation |i,IE0lll3| 

Existing theories of heterogeneous nucleation are phe- 
nomenological and assume that the homogeneous theory 
can be adapted to the heterogeneous case 0, IE S lioj - 
The predictions of these theories have been compared to 
experiment, but the comparisons have been indirect by 
necessity, and the phenomenological nature of the theo- 
ries makes the connection to experiment tenuous. 

The existing theories of heterogeneous nucleation as- 
sume that the homogeneous nucleation process is classi- 
cal, that is, the droplet is assumed to have a distinct vol- 
ume and surface and the structure of the droplet interior 
is the same as the stable phase |E IE • theoretical or 
numerical studies have been undertaken of the effect of 
impurities on nucleation near the pscudospinodal [lllll2l |. 
where in the homogeneous case the droplets are diffuse 
with no sharp surface/volume distinction 11,*^. This 
lack is an important omission because many systems such 
as metals 0| and polymers ^El have a long-range/near- 
mean-field nature with important pscudospinodal effects. 
We will find that heterogeneous nucleation in near-mcan- 
field systems with long-range interactions differs in im- 
portant respects from heterogeneous nucleation in sys- 



tems for which homogeneous nucleation is classical. 

In this letter we discuss our simulations of the effect 
of an isolated impurity and a wall on nucleation. The 
simulations were done on d ~ 2 Ising models at temper- 
ature T — 'iTc/9 with an applied magnetic field of mag- 
nitude h. We used the Metropolis algorithm for about 
100 Monte Carlo steps per spin (mcs) to equilibrate the 
system and then reversed the field. The magnetization 
of the stable phase is negative. As we will discuss, the 
system generally remains in a metastable state for some 
time before the transition to the stable phase. We sim- 
ulated systems with interaction range i? = 1 (nearest 
neighbor) and i? = 10 with linear dimension L = 100 
(unless otherwise stated), and R ~ 20 with L = 200. 
The isolated impurity consists of 5 spins in the shape of 
a -I- sign, fixed in the direction of the stable phase. The 
effect of the impurity on the nucleation rate, the prob- 
ability of heterogeneous and homogeneous nucleation as 
a function of quench depth h, and the structure of the 
critical droplet was studied. 

In Table ^ we show /impurity, the fraction of droplets 
that form on the isolated impurity, and /other, the frac- 
tion of droplets that form on non-impurity sites, even 
though an impurity is present. The number of trials was 
in the range of 100 to 1000. For both i? = 1 and i? = 10 
and shallow quenches, that is, close to the h — coexis- 
tence curve, the droplet forms on the impurity in all of 
our trials. Because the nucleation process is stochastic, 
we expect that a few droplets would nucleate on non- 
impurity sites if we do significantly more trials. As h 
is increased, more droplets form on non-impurity sites. 
This result was found for short-range systems using a 
phase field model 0. 

However, there is a significant qualitative difference for 
deeper quenches between the R = 1 and i? = 10 systems. 
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To see this difference, we plot in Fig.nthe magnetization 
per spin m as a function of time for values of h where 
the fraction of droplet formation events on non-impurity 
spins becomes significant. For R = 10 there is a distinct 
plateau where m remains roughly constant. This plateau 
is associated with the metastable state. However, for 
R = 1 there is no plateau at this quench depth, and thus 
the system is not in a metastable state. 



TABLE I: The fraction of droplet formation events at spins 
away from the impurity, /other, and on the impurity, /impurity, 
for R — 1 and i? = 10 as a function of the quench depth h for 
L = 240. As shown in Fig.0 the droplet formation events for 
deeper quenches correspond to nucleation only for R — 10. 
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FIG. 1: (color online) The magnetization per spin m as a 
function of time (mcs) for R = 1 at h = 0.40 (x) (right axis) 
and R = 10 at h — 1.15 (+) (left axis). Note the flat region 
of m where the i? = 10 system is in metastable equilibrium. 
There is no metastable equilibrium for J? = 1 and h = 0.40. 

In nucleation a critical droplet is a saddle point object, 
which implies that it has an equal probability of growing 
to the stable state or shrinking back to the metastable 
state if the system is perturbed at the time of nucle- 
ation P, Q, 13 ■ We first ran the simulation until a droplet 
formed and the system proceeded to the stable state. The 
spin configurations and the current state of the random 
number generator were saved at various times. We then 
chose an intervention time at which the critical droplet 
might have appeared, made 20 copies of the system, and 
restarted the runs with a different random number seed 



for each copy If the intervention time corresponds 
to the formation of a saddle point object, one half of the 
copies will go to the stable phase at approximately the 
same time and place as in the original run and one half 
will return to the metastable state. If greater than one 
half return to the metastable state, the intervention time 
is too early. If greater than one half proceed to the sta- 
ble phase, the intervention time is in the growth phase. 
Although there will be significant fluctuations with 20 
trials, this number of trials is sufficient for our purpose. 



TABLE II; The number n of interventions (out of 20) for 
which the system goes into the stable phase at approximately 
the same time and place as in the original run; t corresponds 
to the intervention time when the random number seed was 
changed. (The term "no impurity" signifies that no impurity 
was present in the system.) 





t (mcs) 


n 


R — 1 (no impurity) 
h = 0.6 


140 148 151 155 160 


1 5 11 17 20 


R = 1 (with impurity) 
h = 0.4 


70 90 150 250 500 


10 5 15 10 20 


-R = 10 (no impurity) 
h = 1.2 


108 109 110 111 112 


1 8 10 19 20 


R = 10 (with impurity) 
h = 1.16 


198 199 200 201 202 


3 5 15 19 20 



In Table ^ we give the number of interventions that 
proceeded to the stable phase at approximately the same 
time and place as in the original run. Note that a saddle 
point structure is found for the homogeneous and het- 
erogeneous critical droplets in the deeper quenches for 
R — 10. That is, the fraction of interventions that go to 
the stable phase is an increasing function of the time of 
intervention. Similar behavior was found for both values 
of R for shallow quenches, with and without the presence 
of an impurity (data not shown). However, for R = 1 
in the presence of an impurity and quench depths deep 
enough so that we begin to see a non-zero probability of 
droplet formation at non-impurity sites, no saddle point 
structure is found. In this case the fraction of droplets 
that proceed to the stable phase at the same time de- 
creases and then increases as the intervention time in- 
creases, in contrast to the behavior found for i? = 10 and 
for i? = 1 in the absence of an impurity. The lack of a 
saddle point object for i? = 1 and deep quenches is consis- 
tent with the fact that no plateau was found in Fig.^ We 
conclude that the crossover to a significant probability of 
homogeneous nucleation as the quench depth increases is 
found in long-range interaction systems for which the nu- 
cleation is not classical 0. However, in short-range sys- 
tems the "crossover" occurs at quench depths for which 
the system is no longer metastable and the decay process 
cannot be considered to be nucleation llSll. 
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Another quantity of interest is the survival curve, 
which we define as the fraction of systems that remain in 
the metastable state after a given time We perform 

this measurement in a shghtly different but equivalent 
way to that of Heneghan et al. 01 ■ We prepare 1000 sys- 
tems with the same initial condition, but with different 
random number seeds, and then measure the fraction of 
systems, <s>, that are in a metastable state for a given 
value of h after 10^ Monte Carlo steps per spin. 
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FIG. 2: (color online) Survival curves for heterogeneous nu- 
cleation (o), homogeneous nucleation (no impurity present) 
(x), and heterogeneous nucleation on a wall (•) for R = 1. 
The curves for nucleation in the presence of the impurity 
and on the wall are shifted to the right by Ah = 0.174 
and Ah = 0.255, respectively, to make it clearer that their 
shape is qualitatively different from the homogeneous nucle- 
ation curve. 




FIG. 3: (color online) Survival curves for homogeneous nucle- 
ation (x) (no impurity) and heterogeneous nucleation (o) for 
R = 10. The latter curve is shifted to the right by Ah = 0.026 
to make it clearer that the shape of the two curves is almost 
identical. 

In Figs. El and 01 we plot the survival curves for i? = 1 
and R = 10, respectively. Note that for i? = 1 the shape 
of the survival curves for heterogeneous nucleation on 
an impurity and on a wall is qualitatively different from 



homogeneous nucleation (no impurity present) . The wall 
was implemented with periodic boundary conditions in 
one direction and open boundaries in the other, and a 
row of fixed spins at x = 0. In contrast, for R = 10, 
the shape of the survival curves for homogeneous and 
heterogeneous nucleation is very similar (see Fig. |3J|. 

In Fig. ^ the survival curves for R = 10 for homoge- 
neous nucleation (no impurity present) and nucleation on 
a wall are shown. The curves have been shifted to lay on 
top of each other and are almost identical. The results in 
Figs. 13 and pimply that homogeneous nucleation, nucle- 
ation on a wall, and nucleation on an isolated impurity 
have survival curves with similar shapes for long-range 
interactions. This behavior is reminiscent of the exper- 
imental result in Ref. j^l where the survival curves for 
water nucleating on a wall or on a small crystal of Agl 
have similar shapes. The survival curves for R = 20 are 
similar to our results for i? = 10 (see Fig.O. We also did 
simulations without fixed spins at .t = and obtained re- 
sults similar to our other implementation of a wall (data 
not shown). 
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FIG. 4: (color online) Survival curves for homogeneous nu- 
cleation (x) and heterogeneous nucleation (o) on a wall for 
R = 10. The heterogeneous curve is shifted to the right by 
Ah = 0.64. 

For R = 1 and shallow quenches where the nucleation 
is classical, Fletcher showed theoretically that the ra- 
dius of the critical droplet was the same for heteroge- 
neous and homogeneous nucleation. Because the droplet 
interior in classical nucleation is the same as the stable 
phase, the critical droplets have the same structure aside 
from the impurity. However, for R = 10, nucleation takes 
place near the pseudospinodal Hll. [ia|. and the critical 
droplets have a diffuse structure \A Il6l|. The question 
arises as to possible changes in the internal structure due 
to the impurity. In Fig. we plot the density profiles 
of the critical droplets that form on and away from the 
impurity for R = 10. Within the accuracy of our data, 
the structure appears to be the same. 
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FIG. 5: (color online) The survival curves for homogeneous 
(no impurity present) ( ■:) and heterogeneous nucleation (o) 
for R = 20. The heterogeneous curve is shifted to the right 
by Ah = 0.006 to make it clearer that the two curves have 
similar shapes. 



10 



A 
Q. 

V 



0.1 



0.01 




10 

r 



100 



FIG. 6: (color online) The average density <p> of the crit- 
ical droplet as a function of the radius r measured from the 
center of mass for R = 10. (+ signs indicate heterogeneous 
nucleation and • signs indicate homogeneous nucleation.) 



We now summarize our results and discuss their signif- 
icance. First, heterogeneous nucleation droplets appear 
to be saddle point objects, which implies that over a sig- 
nificant range of quench depths, simple modifications of 
homogeneous nucleation theory is possible. Second, the 
common wisdom 0,01 that there is a significant crossover 
from heterogeneous to homogeneous nucleation with in- 
creasing quench depth appears to be partially correct. 
For i? = 10 there is a true crossover because the sys- 
tem is still in the metastable state, whereas for i? = 1 
the "crossover" occurs only after the system has been 
quenched beyond the metastable state. The existence 
of a crossover is important for the theoretical treatment 
of nucleation in systems such as metals where there are 
long-range interactions and the presence of defects 
can be important. 

We found that the survival curves in long-range sys- 
tems have the same shape for both homogeneous and 



heterogeneous nucleation. This behavior is similar to 
the results found experimentally in water |^. However, 
the survival curve shapes for homogeneous and heteroge- 
neous nucleation arc not the same for R = I. 

The connection between the survival curve shape and 
the interaction range is very difficult to make experimen- 
tally because the latter usually cannot be varied system- 
atically. This connection gives us the first experimen- 
tally accessible marker for when a system exhibits nucle- 
ation characteristics associated with long-range interac- 
tions 0,0,^3- That is, if the survival curves have the 
same shape for homogeneous and heterogeneous nucle- 
ation, near- mean-field effects appear to be important. 

We also found that the density profile of the critical 
droplets for heterogeneous and homogeneous nucleation 
is the same within the accuracy of our data for i? = 10 
and that heterogeneous nucleation is a saddle point pro- 
cess in the systems we studied and hence can be described 
by a modified version of homogeneous nucleation. These 
results are significant for the calculation of nucleation 
rates and will be discussed in a future publication. 
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